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(§) Method of developing complex tool shapes. 

@ A method is disclosed for developing the contours of forming tools for fabricating metal members of 
complex shape such as wing skin and fuselage panels for aircraft. A computer simulation is performed of the 
age forming process on a geometrical representation of a member having the material properties of a desired 
metal member. The age forming process includes the steps of: a) overforming an unformed member in a tool 
having a contour of smaller curvature than the contour of the desired member; b) constraining the unformed 
member in the overformed condition; c) applying a thermal aging cycle to the member; d) cooling the 
constrained member following the thermal aging cycle; and e) releasing the constrained member from the 
condition imparted by step (b) and allowing it to spring back to a dimensionally stable condition which defines 
the desired member having a surface contour of complex shape. From the simulation, the contour of smaller 
curvature in step (a) to produce the desired metal member following step (e) is determined. 

FIG. 5 
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The present invention relates generally to a method of developing the contours of fomning tools for 
metal members exhibiting complex shapes and, more particularly, to such a method which utilizes the 
principles of age forming for forming the member being fabricated. The techniques of the present invention 
represent an improvement over those disclosed by H. Brewer and M, Holman in U.S. patent application 
5 Serial No. 07/713.399. filed June 10, 1991, entitled "Method of Tool Development". 

The complex shapes of the contoured members that make up aerospace structures are inherently 
difficult to form. Due to the shapes required by aerodynamics and because of the emphasis on load 
carrying capability combined with weight efficiency, optimized designs are created that require complex 
contours to be produced In high strength metals. Examples of such contoured members would include wing 

70 skin panels, fuselage panels, and structural stiffening elements such as spars and stringers for aircraft 
applications; as well as the shroud, skirt, and tankage members of space launch vehicles. Such members 
are characterized by extreme metal thickness variations and integrally machined features. The criticality of 
design requires precise forming tolerances be maintained without sacrificing the fatigue life, reliability, or 
strength of the member as a result of the forming process chosen. 

75 Conventional forming methods, such as roll forming, brake forming, stretch forming, and peening. are 
cold working processes that achieve permanent deformations through the application of mechanical bending 
and/or stretching. Achieving uniform forming across integrally machined features or abrupt changes in 
thickness may not be possible without specialized tooling or extensive modifications to the forming 
equipment. In some cases, it may not be possible to develop the deformation forces necessary to 

20 accommodate extreme material thicknesses. 

While various machines can handle a wide range of metal thicknesses, it is not practical to form metals 
varying from one extreme of the thickness range to the other, since most machines must be set up prior to 
operating. From this standpoint, skin tapers and recesses that occur within a panel may not be formable. 
Forming applications that have openings or cutouts machined into them may not be formable without 

25 distorting the opening or leaving flat spots in the contour. Other processes are limited by the size of the 
forming machinery and those applications that will fit within the working envelope. Custom equipment for 
larger or smaller applications can be prohibitively costly and inflexible. 

In addition to the physical limitations imposed by part geometry are characteristic traits that result from 
the forming process used. Traits such as strain hardening, residual stresses, and marking accompany many 

30 of the forming processes commonly employed. In some cases these effects can produce desirable 
qualities, such as stress corrosion cracking resistance. Likewise others can produce undesirable qualities, 
such as a negative effect on the fatigue life and reliability of the formed part. The point to be made is that 
each forming process must be carefully matched to the intended application. 

All of the conventional forming processes mentioned have one important disadvantage in common: 

35 each requires the expertise of a skilled operator. With the exception of some processes which have been 
automated to an extent, considerable operator skill is required to obtain tight tolerances; therefore, process 
consistency is low. Part to part variations in contour can result in engineering specified contour rework 
being required on every unit produced. Contour variations that do not require post forming corrections can 
still cause fit-up problems at assembly. Contour variations from part to part create numerous manufacturing 

40 difficulties, each with costly solutions. 

In the recent past, a significant advancement of known techniques for forming complex members while 
maintaining or even improving upon their inherent strength characteristics has been devised. Known as age 
forming, it is a process that offers many solutions to the problems encountered when conventional cold 
forming processes are applied to complex shaped contoured members. During age forming, a part is 

45 restrained to a predetermined tooling contour and precipitation aged. Age forming is a process that utilizes 
the phenomenon of metallurgical stress relaxation during precipitation heat treatment for the purpose of 
converting elastic strain to a plastic state. 

The age forming process may be performed on any of the precipitation heat-treatable, metals and metal 
alloys such as aluminum alloys in the 2xxx. 6xxx, 7xxx, and 8xxx series. 

50 For example, to date, the age forming process of the invention has been successfully employed on at 
least the following aluminum alloys: 
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2xxx Series: 


2014 








2124 




22\A 




2219 




2419 




2090 


6xxx Series: 


6013 




6061 


yxxx Series: 


7075 




7150 




7475 


8xxx Series: 


8090 



Age forming is performed according to standard heat treatment cycles utilized in precipitation hardening 
of alloys. The underlying principles of precipitation heat treating are explained in '•Aluminum Properties and 
Physical Metallurgy". Edited by John E. Hatch, American Society for Metals , Metals Park, Ohio, 1984, pp. 
134-138 and 177-188, which is incorporated herein in its entirety by reference. As a result, suitable 
applications require the final condition of the formed components to be in an artificially aged temper. Every 
end use of a structure must be reviewed in light of the property changes that occur as a result of artificial 
aging. In some cases, the mechanical properties associated with an artificially aged temper may not be 
suitable for an Intended application. As an example, aluminum alloy 2024 loses fracture toughness as it is 
artificially aged from the T3 to the T8 temper. This change presents a barrier to age forming applications 
where fracture toughness is a key design element, such as lower wing skins and fuselage panels for 
aircraft. Material and/or design changes are required in these cases to allow for the utilization of age 
forming. In other cases, age forming allows the added benefit of being able to produce contours in a 
strengthened temper, without developing high levels of residual stress within the component. An example of 
this feature is provided when aluminum alloy 7150 is age formed from the soft W temper to the hardened 
T6 temper. 

More recently, the conventional age forming process has been modified and substantially improved 
through the use of the autoclave. The autoclave is a computer controlled pressure vessel, with the added 
benefit of being a certifiable source for heat treating aluminum. Age forming has traditionally been 
performed in a furnace, where a mechanical means of constraining the part to the predetermined forming 
shape is required. The autoclave offers the advantage of using vacuum and internal pressure to obtain the 
desired contour. Since pressure acts uniformly about the surface of the part, integrally machined features 
receive the same deformation force as the rest of the panel. Another important advantage is that the 
forming pressure is distributed about the entire surface area of the part. Therefore, a small differential 
pressure can equate to many tons of applied force when acting over a large surface. Most conventional 
processes concentrate the forming forces over a small area, thereby restricting the total available 
deformation force. 

The autoclave is computer controlled allowing high levels of process consistency and accuracy. 
Computer control allows the process to be operator independent. A separate computerized system closely 
monitors and records the pressure and temperature within the autoclave providing traceability and process 
verification. These two features inherently endow autoclave age forming with high levels of process 
consistency and accuracy. Each panel receives the same processing; consequently, repeatability is 
ensured. It is this feature that makes the process adjustable. The tooting contour is "fine tuned" until the 
desired results are obtained. 

Tooling for the autoclave is designed according to the springback anticipated for the application. 
Springback refers to the tendency for a member being formed to return to some shape intermediate its 
original shape and that of the tool to which it is subjected during heat treatment. This phenomenon will be 
discussed at length below. Forming tools are designed with removable contour boards and other features 
that allow for rapid contour modifications. Unlike other forming processes, age forming does not typically 
allow for multiple forming iterations to be performed upon the same piece. Age forming is a heat treatment 
process; therefore, running a part more than once could result in over aging the material. Until the tooling 
contour is finalized, contour corrections must be performed by another forming process. Once the final tool 
contour is reached, secondary corrective forming processes are not necessary. 
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This inability to repeat the heat treatment process on a member being fabricated requires that it be 
scrapped if it exhibits an inconrect final contour and the procedure repeated with a new member. The cost 
of labor and materials for such necessarily repeated iterations of the process have led to the methods of the 
present invention. 

5 According to the invention, there is provided a method of forming a desired metal member having a 
surface contour of complex shape from an unformed member comprising the steps of: 

a) overforming the unformed member In a forming tool having a contour of smaller curvature than the 
contour of the desired member; 

b) constraining the unformed member in the overformed condition; 
70 c) applying a thermal aging cycle to the member; 

d) cooling the constrained member following the thermal aging cycle; 

c) releasing the constrained member from the condition imparted by step (b) and allowing it to spring 
back to a dimensionally stable condition which defines the desired member having a surface contour of 
complex shape; 

75 (f) simulating the application of steps (a) through (e) to a geometrical representation of a member having 
the material properties of the desired metal member; and 

(g) determining from step (f) the contour of smaller curvature in step (a) to produce the desired metal 
member following step (e). 

Thus, a method is disclosed for developing the contour of tools employed for forming metal members 
20 exhibiting complex shapes. The members are precipitation, heat-treatable, metals or metal alloys which are 
autoclave age formed. The resulting member is formed to the desired contour and, simultaneously, is heat 
treated to reduce residual stresses while improving its strength characteristics. The invention Is particularly 
concerned with a new tooling contour prediction method which is based upon the application of a unified 
viscoplastic model simulating the age forming process. 
25 The method of the invention assures proper results on the first occasion the tool is used, thereby 
resulting in considerable savings of labor and material. 

There now follows a description of a preferred embodiment of the invention and of preferred method 
steps in accordance with the invention, with reference being made to the accompanying drawings, in which: 
Fig. 1 is a diagrammatic side elevation view illustrative of stress distribution in a constant thickness bar 
30 being subjected to pure bending for purposes of explanation of the invention; 

Rg. 2 is a stress-strain graph illustrating the relationship between stress and strain in the outermost layer 
of material of the bar of Rg. 1 during a cold mechanical forming process, depicting both the elastic 
range of the material and the deformation in the material after It has been stressed beyond the yield 
strength of the material; 

35 Rg. 3A Illustrates a stress-strain graph, similar to Rg. 2. but indicating the result of an age forming 
process performed within the elastic range of the material; 

Rg. 3B is a stress-strain graph, similar to Rg. 2, but indicating the result of an age forming process 
performed when initial loading exceeds the yield point of the material; 

Rg. 4 Is a perspective view, exploded, illustrating tooling for autoclave age forming a member such as 
40 the bar of Rg. 1; 

Rg. 5 is a detail cross section view illustrating the items shown in Rg. 4 within an autoclave; 

Rgs. 6A, 6B, 6C are successive diagrammatic detail end elevation views, partially in section, illustrating 

successive steps of the age forming method of the invention; 

Rg. 7 is a basic flow chart of the simulation model of the invention; 
45 Rg. 8 is a diagrammatic cross section view of a forming tool and unstressed part thereon presented on 

Cartesian coordinates; 

Rg. 9 is a graph depicting the predicted stress distributlonon the central cross section of a member for a 
range of pressures to exhibit the slopes characteristics of elastic plastic response; 
Rg. 10A diagrammatically represents a tool reference plane; 
50 Rg. 10B diagrammatically represents a desired metal member cut by a plurality of parallel spaced apart 
planes perpendicular to the tool reference plane of Fig. IDA; 

Rg. 11 is a diagrammatic representation. In section, of a cross sectional element resulting from the 
interface of one of the perpendicular planes in Rg. 108 and the desired metal member; 
Rg. 1 1 B is a diagrammatic representation of a plurality of segments resulting from the cross sectional 
55 element of Rg. 11 A; 

Rg. 12 is a process flow chart presenting the method of the invention in more specific terms; 

Rg. 13 presents a graphic illustration of the method by which a smooth continuous curve is achieved 

utilizing the present invention; 
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Fig. 14 is a detail perspective view of a part of a tool resulting from the method of the present invention; 
Fig. 15 is a diagrammatic view of a flattened cross sectional element enabling its measurement on 
Cartesian coordinates; 

Fig. 16 is a diagrammatic representation depicting a comparison between the shape of a simulated age 
5 formed member and the shape of a cross sectional element; and 

Fig. 17 is a graph depicting, for a range of thicknesses of a member to be formed, the relationship 
between the radius of curvature of a forming tool and the radius of curvature of the formed part. 
In order to gain a better understanding of the phenomena behind the age forming process of the 
invention, it is well to separately consider and analyze the forming mechanisms at work during the age 
70 forming process. This effort can begin by analyzing mechanical forming versus age forming in terms of 
stress distribution found within the cross section of a specimen undergoing forming. Another tool desirably 
utilized for analysis Is a stress-strain curve representing the outside layer of fibers of a specimen 
undergoing forming. Through the use of these tools, a clearer picture can be obtained as to how each 
forming method works to form a piece of material. 
75 Considering the stress distribution throughout a part 20. depicted for simplicity in Fig. 1 as a constant 
thickness bar of rectangular cross section, allows a comparison to be drawn between different forming 
mechanisms. As a force F is applied to the simply supported bar to cause it to assume a radius, stresses 
diagrammatically Indicated at 22 are distributed throughout the thickness of the bar. A neutral surface 24 
experiences no stress due to pure bending while the outside fibers experience the greatest stress. A 
20 concave side 26 of the bar experiences compressive stresses while a convex side 28 of the bar 
experiences tensile stresses. According to Hooke's Law, stress is directly proportional to the strain that is 
experienced when it is within the elastic range of the material. The proportionality constant is known as the 
modulus of elasticity and is dependent upon material and temperature. The strain experienced by the fibers 
across the thickness of a specimen depends upon the distance of a particular layer of fibers from the 
2s neutral surface. 

If the stress induced throughout the bar stays within the elastic range of the material, the bar will return 
to Its original flat configuration with no forming taking place once it is released. Therefore, if the bar is to 
retain a contour and be formed without the aid of thermal stress relaxation, a significant amount of fibers 
within the material must be stressed beyond their yield point. The stress-strain curve 30 In Fig. 2 can be 

30 used to examine the action involved in forming. The case of imparting a radius to a flat bar shaped part is 
not strictly a tensile application; rather it is one of bending. Therefore, in reality, the use of a stress-strain 
curve is only applicable to a single layer of material at a given distance from the neutral surface. 
Nevertheless, it serves the purpose of Illustrating the differences between cold mechanical forming and age 
forming. For example, the stress-strain curve 30 In Fig. 2 illustrates cold mechanical forming of the bar 20 

35 of Fig. 1 subjected to bending stresses. 

Consider the outermost layer of material on what will become the convex side 28 of the bar. Initially the 
bar Is flat and in a stress free state. As the bar is reconfigured to assume a radius, the fibers in the outside 
surface layer are strained which induces stress. This is illustrated by a stress distribution line 32 (Fig. 1) and 
by the stress-strain curve (Fig. 2) beginning at the origin. The linear portion of the curve, which defines the 

40 modulus of elasticity, or Young's modulus, for the particular alloy of the bar 20. continues until the stress 
level reaches the yield strength 36 of the material. If the bar is released at any point prior to inducing a 
stress greater than the yield strength 36, it will unload along this same line and return to a flat (i.e., strain 
free) condition. Once a layer of material Is stressed beyond its yield point, the relationship between stress 
and strain Is no longer directly proportional (i.e., it is no longer linear). If at this point the bar Is released, it 

45 will unload along a line 38 that has the same slope as the linear portion 34 of the load curve 30 but will be 
offset from the original load line 34 Indicating a retained strain 40. The slope is equal to the modulus of 
elasticity as previously noted. The resulting retained strain 40, referred to as plastic strain, indicates that 
permanent deformation has taken place. 

Age forming forms a structure by taking advantage of the stress relaxation phenomena associated with 

50 artificial aging. The age forming concept Is Illustrated by the stress-strain curves In Figs. 3A and 3B. Fig. 3A 
depicts a specimen initially stressed below the material's yield strength and Fig. 3B depicts a specimen 
initially stressed beyond the material's yield strength. Again, consider the outside layer of fibers on what will 
become the convex side of a formed member, such as convex side 28 of the bar 20 of Fig. 1 . These fibers 
will experience tensile stresses. As the member Is strained as indicated by a line 42 (Fig. 3A), the stress 

55 level Increases proportionally. The vertical distance oa (Fig. 3A) represents the amount of stress exper- 
ienced by the fibers of the member while the horizontal distance €i represents the amount of strain 
experienced. Upon reaching a particular radius, the member is held at this constant strain level (as at point 
44) and the artificial aging cycle is applied. Due to the metallurgical stress relaxation resulting from the 
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materials* exposure to temperature, the stress level reduces even though the strain remains constant. The 
amount of stress relaxation that occurs, as indicated at Ob, depends upon the material and its related aging 
temperature as well as the initial level of stress induced. The rate of stress relaxation is greatly enhanced by 
a higher Initial stress level and by a higher aging temperature. However, these factors are limited by the 

5 temperature permitted by the selected aging cycle. Once the aging is complete, the member is cooled and 
released from Its constraints. This allows the member to spring back and physically relax the remaining 
induced stress. The vertical distance oc (Fig. 3A) represents the amount of stress relaxed during spring 
back while the horizontal distance fa represents the change in strain. Since strain changes, the shape of the 
specimen also changes. In this case the specimen is held In contact with the smaller radius of a forming 

10 tool and. upon release and following spring back, assumes a larger radius. An amount of strain €2 is retained 
by the member indicating permanent deformation. 

In Rg. 3A. the practice of age forming has been illustrated within the elastic range of the material. It is 
in this region that the distinction between age forming and cold mechanical forming is most evident. 
However, the same principles apply within the inelastic range (above yield) as depicted In Fig. 3B. The most 

75 notable difference between age forming a specimen stressed within the elastic range versus one stressed 
within the inelastic range is best viewed by considering the action along the strain (horizontal) axis. In a 
specimen stressed to within the inelastic range, the retained strain €15 (Fig. 3B) is composed of two 
components. A portion of the retained strain C12 results simply due to stressing the specimen beyond the 
yield point of the material- In Fig. 3B. paint xx represents the specimen initially reconfigured to the shape of 

20 the forming tool prior to exposure to the aging cycle. At this point, the level of stress Is beyond the yield 
strength of the material. The yield strength is illustrated on the stress-strain curve 42A by point yy. If the 
specimen being formed were to be released at point xx. prior to exposure to the elevated temperature of 
the aging cycle, some retained strain €12 would be exhibited simply because a portion of the material has 
yielded. This is unlike the specimen illustrated in Fig. 3A in which the specimen would return to a flat 

25 unstrained condition if released prior to elevated temperature exposure. The total retained strain €15 of Fig. 
3B. therefore, is a combination of the retained strain €12 due to yielding of the material and the retained 
strain cia due to metallurgical stress relaxation. 

In either the elastic or inelastic range, age forming allows permanent deformation to be achieved with 
lower levels of applied stress than cold mechanical forming. Because of the way that cold mechanical 

30 forming works, residual stress levels within formed parts can be quite high. It is here that age forming 
presents significant advantages. Rrst, the applied stress level required for forming is lower; and secondly, 
stress relaxation occurs during aging, lowering it even more while the part is held at a constant strain. After 
release from the forming tool, the age formed part relaxes the remaining Induced stress, which is 
significantly lower than it was at the start of the aging cycle. The result is that the age formed part has the 
35 same permanent deformation as the mechanically formed part, but with much lower levels of residual 
stress. 

The amount of stress relaxation experienced by a member during forming becomes the key to 
determining the amount of springback the member will experience following age forming. Predicting 
springback is the fundamental requirement to taking advantage of the age forming method. Knowledge of 

40 springback is needed to accurately determine forming tool contours. 

For a brief initial explanation of the autoclave age forming process utilized for purposes of the invention, 
turn now to Figs. 4 and 5. An autoclave 50 (Fig. 5) includes a generally thick-walled cylindrical vessel 52 
which may typically be capable of withstanding pressures up to 200 psi. total vacuum, and temperatures up 
to 600**F. With this apparatus, as diagrammatically seen in Fig. 6. the part 20 is forced from an initial 

45 unformed condition (Fig. 6A) into intimate contact with the contoured surface 53 of a concave die 54 (Fig. 
6B) receivable in a cavity 56 of an autoclave forming tool 58. This is accomplished by covering the top of 
the part 20, die 54. and forming tool cavity 56 with a temperature resistant vacuum blanket 60. sealing the 
edges of the blanket, drawing a vacuum through a plurality of vacuum ports 62 (Fig. 4) on the tool cavity 
beneath the part, and, if desired, also applying pressure to the upper surface of the part. A sealing frame 64 

50 is removably mounted on the forming tool 58 to maintain the positioning of the vacuum blanket 60. The 
vacuum pulled underneath the part ensures that trapped air will not prevent it from obtaining total contact 
with the forming tool. The forming tool contour is designed to overform the part, allowing for springback. As 
noted above, pressure may be optionally applied to the part as indicated by arrows 66 to assure firm and 
continuous coextensive engagement of the die 54 by the part 20. 

55 Up to this point, temperature has not been applied to the part, so that unless the bending stress applied 
has exceeded the yield point of the material, no permanent deformation has been achieved and the part is 
still within the elastic range of the stress strain diagram. This condition provides the most significant feature 
of age forming, since it can be performed at lower applied stress levels than conventional forming 
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techniques. If the part were released from the vacuum and pressure holding it to the tool, it would 
essentially springback to its initial flat condition (Fig. 6A). However, with the application of heat at 
appropriate temperatures for appropriate periods of time, the part will, after the forming process and after its 
release from the tool, spring back to an intermediate position as indicated in Fig. 6C. 
5 The foregoing presents an early construction of an autoclave tool suitable for the process of the 
invention. However, it is not all inclusive. More recently, tools have been constructed with a skeleton 
framework of contoured boards covered by a contoured aluminum skin or caul plate. The pressure 
differential is created between the top of the panel and the caul sheet. The contour boards are not exposed 
to the pressure differential, except for those forces transmitted through the caul. A sealing frame is no 

10 longer employed to seal the vacuum bag to the tool. Instead, the vacuum seal is now maintained by 
adhesively attaching the bag to the surface of the caul with a temperature resistant putty. The newer tooling 
is simple, light-weight, and less costly to build. Nor does the tooling have to be concave; it can just as 
easily be convex. Also, production tools are not generally cylindrical, although individual contours are 
constructed of circular segments. While vacuum and pressure are preferably employed to obtain the 

15 appropriate applied strain, purely mechanical expedients, such as matched dies or clamps, may also be 
used. Much of the tooling is simply a function of the desire to use a pressure differential for forming. Age 
forming itself can be employed in both autoclaves and furnaces using both pressure and mechanical 
means. The method for developing the forming tool contour is the same, regardless of whether a 
pressurized autoclave tool or a mechanically clamped furnace tool is desired. Springback is calculated as a 

20 function of the material, its thickness, and the final contour desired only. Regardless of whether age forming 
is performed in a furnace or autoclave, the material's response to aging remains the same. 

Until the advent of the present invention, springback was defined as the difference between the chord 
height of the tool and the chord height of the formed specimen. However, It was found that this method was 
very restrictive and limited to predicting the springback of a constant thickness bar specimen formed to a 

25 radius. The old method was based purely on the percent change in chord height. The stress-strain curve 
was not used. This method was improved by using the stress relaxation curve and strain retention curve 
prediction method as indicated in U.S. Application Serial No. 07/713.399 recited above, the disclosure of 
which is hereby incorporated herein in Its entirety by reference. However, the improved method, just noted, 
is based on experimental observations and was limited to the range of test data that was used. 

30 A new springback prediction method which forms the basis of the present invention is based upon the 
application of a unified viscoplastic model to simulate the age forming process, providing a much more 
complete analytical device than previously available to the tool designer. The age forming method can be 
broken down into its various stages: loading, stress relaxation, and springback. A basic flow chart depicting 
this method is presented In Fig. 7. In the present method, equations representing the condition model are 

35 used to more fully describe what is physically and metallurgically happening to the material being formed. 
These equations attempt to describe the laws governing the physical nature of the material and changes 
taking place during the age forming process. 

They represent physical phenomena such as: elastic strain, inelastic strain, stress relaxation, creep, and 
the like, and the history of time dependent load application and temperature exposure. Unique constants are 

40 required to accurately represent specific materials. The constants are determined by manipulating the 
constitutive equations until they represent the age forming process physically observed in test specimens. 

Once determined, the constants in conjunction with the constitutive model fully represent the material at 
hand as it undergoes the age forming process. 

Theoretically, any model geometry can then be analyzed to determine needed age forming tool 

45 contours. More properly, this method is a modelling and simulation technique rather than a prediction 
technique. 

Mathematical modelling and simulation of age forming is flexible and Incorporates material properties 
and part geometry in an appropriate format. The model uses this information to obtain the desired contour 
of a specimen being formed and to predict the residual stress in that specimen. Integrating materials, as 

so represented by the constitutive model, and geometry into the model for the forming method allows it to be 
adaptable to different combinations of part configuration and metal alloy. 

The benefits of a mathematical modelling and simulation of the age fonming method relate to the ability 
to know the degree of deformation required to compensate for material springback and the characteristic 
forming tendencies associated with a specific part configuration. The main benefit would be to analytically 

55 determine the forming parameters thereby eliminating the need for developing costly and time consuming 
empirical data. 

The procedure to obtain the basis for the present invention will be divided into the following sections: 
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INTRODUCTION 

In autoclave age forming, a metal part is heated to a temperature sufficient for it to exhibit creep and 
stress relaxation. It is then subjected to pressure which forces it against a tool or mold. Creep may facilitate 
5 this process. Once it contacts the tool, the part is held in place by the pressure for a period of time to allow 
stress relaxation to relieve the stresses produced by forming. Finally the part is released and it partially 
springs back to a shape somewhere between its original shape and the tool shape. 

From the above description, it can be seen that the process is quite complex. A variety of mechanical 
phenomena including elastic deformation, plastic deformation, creep, and stress relaxation are potentially 
10 involved. This complexity of mechanical response makes the use of a unified model of elastic and inelastic 
behavior attractive. 

CONSTITUTIVE MODEL 

15 A mechanical model capable of a unified description of yielding, creep, and stress relaxation is needed. 
The constitutive model also should have structure involving evolution equations for measures of inelastic 
strain and a set of intemal variables. Microscopic concepts can be used in determining the forms of the 
equations governing the evolution of the internal variables. The evolutionary nature of the method of the 
invention allows a smooth transition from elastic to Inelastic behavior and thus eliminates the need for the 

20 explicit inclusion of a yield criterion. After Investigating several possibilities, the Miller and Sherby 
constitutive model (see Miller, A.K., and Sherby, O.D.. "A Simplified Phenomenological Model for Non- 
Elastic Deformation: Prediction of Pure Aluminum Behavior and Incorporation of Solute Strengthening 
Effects," Acta Metallurglca . Vol. 26, 1978, pp. 289-304) was finally selected. The chosen model Is 
potentially well suited for the simulation of a complicated process such as age forming in which the relative 

25 contributions of such phenomena as yielding, creep, and stress relaxation are initially unknown. In a beam 
model, only the one dimensional forms of the Miller and Sherby equations are needed. These are In a form 
and notation somewhat different from that used in Miller and Sherby. for convenience, but not substantially 
different: 

30 a + (E/T)(sinh(a2/(E2f)))" = Ef (1) 

f = (a/r)((sinh(a2/(E2f)))"-(sinh()820)") (2) 
f(0) = fo (3) 

35 

where a superposed dot indicates differentiation with respect to time (t). 

The symbols a. and f represent, respectively, stress, total strain, and an internal variable which is a 
dimensionless form of the quantity called the drag stress in Miller and Sherby. The inclusion of the internal 
variable f provides a mechanism for roughly modelling the internal structural changes responsible for 
40 inelastic behavior. The quantities a, /8. fo and n are material constants while E and t are temperature 
dependent material properties. The physical interpretations of a, fi, fo. and n are discussed by Miller and 
Sherby. The quantity E Is the familiar Young's modulus of elasticity, and the quantity t Is a relaxation time. 

In order to utilize the Miller and Sherby model to simulate the behavior of a specific material, a set of 
material constants appropriate for that material must be determined. In the present work this was done with 
45 reference to uniaxial constant strain rate tension tests, uniaxial constant stress creep tests, and uniaxial 
constant strain stress relaxation tests which are conducted for the desired alloy. 

All the tests mentioned above Involve uniaxial loading under quaslstatic conditions. For these cir- 
cumstances, the quantities o, and f can be regarded as constant throughout the specimen. The equations 
(1) and (2) are a set of two simultaneous differential equations which can be solved to determine either a(t) 
50 and f(t) when €(t) is given (as in the tension and stress relaxation tests) or ^(t) and f(t) when o(t) is given (as 
in the creep tests). The problem of carrying out the solution of equations (1) and (2) numerically is 
complicated by the numerical stiffness of the system. 

In the present work it was decided to employ the Euler two point forward difference method and to deal 
with numerical stiffness by using a very small constant step size. This procedure was verified by 
55 reproducing several simulations reported by Miller and Sherby for pure aluminum using the material 
constants given in the first column of Table 1. 
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Table 1 



Material Constants 




Pure Aluminum 


7075 Aluminum Alloy 


a 


4.70x10' 


5.00 X 103 




5.50 X 10^ 


3.75 X 102 


n 


4.70 


5.00 


fo 


3.30 X 10-3 


3.00 X 10-3 


r 


2.97 X 10'* sec 


6.06 X 10' sec 


E 


6.96 X 10' kpa 


6.55 X 10' kpa 



It was found that any sudden application of stress or strain produced numerical difficulties. Thus creep 
tests were simulated by increasing the stress to its final constant value over a time which was small 
compared to the duration of the test and stress relaxation tests were simulated by increasing the strain to its 
final constant value over a time that was small compared to the duration of the test. 

Because of the complicated structure of equations (1) and (2), it is impossible to identify each of the 
material coefficients with a unique aspect of the material and response. Thus, the determination of the 
coefficients involves a trial and error process based on a large number of simulations. In the present work a 
procedure similar to that described by Miller and Sherby was used to obtain initial estimates of the 
constants. This procedure is based on data from uniaxial creep tests and uniaxial tension tests. These 
estimates were then refined using the available data for uniaxial stress relaxation tests. For example, the 
values ultimately selected as appropriate for age forming simulations of 7075 aluminum alloy are given In 
the second column of Table 1. 

BEAM MODEL 

A part could be represented by a simply supported beam subjected to uniform pressure. It is assumed 
that the contact between the tool and the part is established everywhere simultaneously and that the contact 
pressure can be treated as uniform. 

It should be noted that, consistent with the above discussion, during the contact phase the pressure 
appearing in the subsequent equations represents the difference between the applied pressure and the 
contact pressure. The beam model will first be developed in a form suitable for the solution of any bending 
problem and then specialized to the age forming problem. 

To develop the appropriate equations for quasistatic bending let the x,y plane t>e the plane of bending 
with the X axis being horizontal and the y axis vertically upward (see Fig. 8). The undeformed neutral axis 
corres|X)nds to y = 0. The cross section of the beam is characterized by an area A and a moment of inertia 
about the neutral axis 



(4) 

The bending moment Is defined to be 



M(x,t) = - I I aydA 

(5) 

The vertical displacement of the neutral surface is v(x,t). 

In the beam model the stresses and strains can vary with x and y as well as t. Thus, the superposed 
dots appearing in equations (1) and (2) must be replaced by partial derivatives with respect to t. Doing this 
and substituting the geometric relationship 



9 



EP 0 527 570 A1 



€ = -y/R(x.t) = -yjr(x.t) (6) 

(R being the radius of curvature of the neutral surface and x being its curvature, that is. 1/R) into equations 
5 (1) and (2) yields 

d,a + (E/r)(sinh(o2/(E2f)))" = -Eyd,* (7) 

d,f = (a/T)((sinh(a2/(E2f)))"-(sinh(/52f))") (8) 



10 



15 



25 



30 



Integrating the product of y with equation (7) over the cross section and using equations (4) and (5) yields 
the additional equation 



^iM - (E/t) j j (sinh(aW0)fydA=EiaiK 



(9) 

20 For small deflections, the relationship between curvature and displacement is 



2 

K = axv (10) 



For statically determinate problems, equations (7-10) are sufficient to define the problem. For statically 
indeterminate problems, they must be supplemented by the quasistatic equilibrium equation 

2 

q = axM 

35 where q is the upward loading per unit length. It is to be observed that for t = » one integration of equation 
(9) with respect to t produces the usual moment curvature equation for elastic being jr = M/(EI) and 
substituting this into one integration of equation (7) with respect to t yields the usual elastic flexure stress 
distribution a = -My/I. 

The integral appearing in equation (9) illustrates the fact that, without further assumptions, a beam (or 
40 plate or shell) theory cannot be rendered entirely independent of the thickness coordinate when a unified 
model of elastic/inelastic response is employed. In the present situation this is because the variations of o 
and f with y are not given but must be determined as part of the solution. 

The solution of equations (7-11) is a formidable problem. In the age forming experiments being 
simulated it was found that the formed parts closely resembled circular arcs. It was decided, therefore, to 
45 simplify the model by assuming that equation (10) could be replaced with the appropriate equation based 
on the geometry of a circular arc. Doing this allows the solution to be carried out at each cross section 
independently. In the present work only the central cross section, the location of the maximum stresses and 
displacements, was considered. 

The age forming model involves a statically determinate simply supported rectangular beam (of length 
50 t. width b. and thickness h) loaded by a uniform downward pressure p (note: q = -pb). Thus the central 
moment is 

M = pbP/8 (12) 

55 and the moment of inertia is 

I = bh3/12 (13) 
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Suppressing the dependence on x in equations (7-9) (since x is fixed at 1/2) and using equations (12) and 
(13) one obtains 

6to + (E/r(sinh(a2/E20))" = -Eyi (1 4) 

d,f= (a/r)((sinh(o2/E2f)))"-(sinh(/82f))") (15) 



Mm 

FK16E/(Tfi)J (8inh(aW0))Vy = 2Eh^ii/(3ft 



(16) 



Equations (14-16) are a set of differential equations which can be used to determine o(y,t) and f(y,t) on the 
central cross section together with either x(t) (if p(t) is given) or p(t) (if x(t) is given). The equation 

V = R(1-(1-(t(2R))2)^'2). ^ 

(see Fig. 8) is used to determine the central deflection. 

Before attempting to simulate the entire age forming process, simulations are made off the process of 
deforming a beam by a pressure which starts at zero and increases linearly with time. The material 
properties for pure aluminum given previously are employed and the dimensions are assumed to be I 
= 76.2cm (30 in.), b = 7.62cm (3 In.), h = 1.27cm (0.5 In.). 

The simulations are carried out as follows. In the upper half of the cross section. N points are chosen 

and the quantities oi, ci. and fi. I = 1 ,2 N are associated with these points. The initial state is assumed to 

be quiescent, that is 

a, = c, = O.f, = fo;i = 1.2 N. (18) 

Since p(t) was known in this case, the governing equations were arranged in the form 

Oi = -((E/T)(sinh(di/(E^fi))f + 3^p/(2h^) 
-(24E/ch^)| (8inh(aWf))fydy)(y/h)); i = 1,2,...,N 

(19) 

fi = (a/x)((sinh(o?/(E?fi)))".(8inh(P^fi))"); i a 1,2,...N (20) 



K = (3lV(2Eh^) - (24/(xh^) (sinh(aV0))"ydy)/h 

h (21) 



and these 2N + 1 ordinary differential equations are integrated forward by the well known Euler method to 
determine ai(t), fi(t), and x(t). It Is interesting to note that, since the Euler method is being used, the integrals 
appearing in equations (19) and (21) are always evaluated at the previous time. Thus, no iteration is 
required. This may not be true for many, more complex, numerical methods. 

The capability of the present beam model is illustrated in Fig. 9 which shows the predicted stress 
distributions on the central cross section. For the lowest three pressures the stress distributions are linear 
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which is characteristic of elastic response. Beginning with the fourth pressure the stress profiles exhibit the 
slopes characteristic of elastic plastic response. The yield stress is approximately 10 Mpa which is 
consistent with data reported by Miller and Sherby. The yielded region is seen to occupy an increasingly 
large portion of the cross section as the applied pressure is increased. For the largest pressure shown, 
virtually the entire cross section has yielded. It is felt that this simulation of the transition from elastic to 
plastic behavior (yielding) is quite good. 

AGE FORMING SIMULATIONS 

The age forming process consists of three stages as discussed earlier. These will be refenred to herein 
as loading, stress relaxation, and springback. 

The loading stage is simulated exactly as described at the end of the previous section. It is terminated 
when the part is predicted to be in contact with the tool, that is. when R = 1/x is found to be equal to the tool 
radius Rf. 

To simulate the stress relaxation stage the curvature of the part is held constant at the tool curvature 
(thus making x =0) and equations (19-21) are rearranged to read 



Equations (22) and (23) consist of N sets of two coupled equations. These are integrated fon^vard by the 
Euler method to determine oi(t) and f|(t). Then equation (24) is integrated forward by the Euler method to 
find p(t). The stress relaxation stage is terminated after 24 hours. 

It is to be recalled that in the present idealized model the quantity p represents the difference between 
the forming pressure and the distributed reaction between the part and the tool. This quantity decreases as 
the stress relaxation stage proceeds. At the end of the stress relaxation stage the forming pressure is 
removed. This relieves the contact between part and tool, thus causing the distributed reaction between part 
and tool to vanish. In the simulations, therefore, the springback stage is characterized by reducing p to zero 
at a constant rate and then holding it there. The simulations are carried out using equations (19-21) just as 
for the loading stage. It has been found that chances in x virtually cease once p reaches zero (that is R = 1/x 
at the end of the unloading stage is equal to the part radius. Rp). The corresponding value of R = 1/x is used 
to characterize the amount of springback. It is also to be recalled that in the actual experiments the part is 
cooled before the pressure Is removed. In the simulations, this cooling is accounted for by the appropriate 
changes in E and r in the unloading stage. 

Determining the proper amount of overform a metal member must be subjected to during age forming 
is critical to successfully forming the desired metal member. The more complex the desired member, 
whether in contour or in thickness variation, the more complex the forming tool contour must be. The unified 
viscoplastic model which has been described simulates the behavior of a metal item when subjected to age 
forming. Therefore, the model can become a key component to an overall scheme for predicting an age 
forming tool's contour. The following presents three schemes for such predictions. In each scheme, a 
complex shaped metal member of desired contour is analyzed and information from that analysis is used in 
conjunction with the simulation model to produce a forming tool contour. 

For example, consider a wing skin panel as the metal member requiring form. Thickness varies 
throughout the panel as required by its intended purpose as a wing covering. Also the shape to which it 



6i = (E/x)(sinh(o?/(Efi)))"; i = 1,2,,..,N 



(22) 



fi = (a/t)((sinh(a?/(E^fi)))" - (sinh(|J fi))"); i = 1.2,...,N 



(23) 




(24) 
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must conform is complex as defined by the airfoil shape. In order to predict the required forming tool 
contour, the shape of the desired metal member must be known. A scheme for analyzing the shape and 
thickness variations of the desired member has been developed which simulates subjecting the member to 
a plurality of two-dimensional cuts or slices. First, a geometrical representation of the desired member is 

5 created in an electronic three-dimensional graphics package such as that sold under the trademark 
"CATIA". Here the shape of the desired member can be fully described in terms of an x. y, z coordinate 
system. Within this medium, a plane of reference is established. This reference plane is labelled as the Tool 
Reference Plane (TRP) 201 (Fig. 10A). This plane represents a desired or contoured metal member 203 in 
a flat condition prior to forming. The desired metal member in its final configuration Is illustrated in Fig. lOB. 

10 With TRP 201 established, conveniently spaced apart planes 205, 207. 209 perpendicular to the TRP 201 
are used to cut or slice the desired metal member across its contour By means of these planes 205, 207, 
209. a desired metal member 203 can be defined in a two-dimensional and cross sectional format as seen 
in Rg. 11 A. Each scheme for determining the forming tool contour utilizes the resulting cross sectional 
elements 211. 

15 Each cross section as represented by its associated cross sectional element 211 provides information 
concerning shape and thickness variation. For the first scheme of determining the forming tool contour, 
each cross sectional element is divided Into segments 213, 215. 217. 219, respectively, containing 
substantially similar thickness (hi.h2,h3.h4. respectively) and substantially similar radius (Ri, Ra, R3. R*. 
respectively) combinations, each with a corresponding length (h, la, I3. U. respectively) as shown in Fig. 

20 11 B. Depending upon the cross section, several segments may come from a single cross sectional element. 
The next step is to determine what forming tool radius would be required to produce the shape of the 
segments descrit>ed above. 

At this point the simulation model is used. The simulation model simulates age forming as applied to a 
flat specimen being formed upon a specifically shaped tool. The model accurately predicts the final shape 

26 of the specimen having undergone the age forming method. In order to fully utilize the simulation, it is 
necessary to vary the shape of the forming tool and compare the results of the predicted final shape with 
the shape of the segment being sought. For each segment, the desired radius sought in the segment, its 
length, and its thickness are supplied to the simulation model The model creates a member having a length 
and a thickness equivalent to the segment being analyzed. An initial tool radius which is smaller than the 

30 radius of the segment being analyzed is supplied to the simulation model. The model then simulates the 
age forming of the created member as if formed against this initial tool radius. The final radius produced by 
the simulation is then compared to the desired radius of the segment being analyzed. Predetermined values 
are established in the simulation model as being acceptable differences between simulated and desired 
values. If the differences between the desired and simulated radii are not acceptable, then the simulated 

35 forming tool radius is adjusted according to these differences. The simulation is performed using the new 
tool radius, and comparisons are made again. The scenario continues in an iterative fashion until the 
differences between the simulated and desired radii are within acceptable levels. Once the differences are 
acceptable, the forming tool radius utilized by the simulation is outputed as the forming tool radius required 
to age form a member of the specified thickness to the desired final radius. This simulation is continued for 

40 each different thickness and radius combination represented by each sement. A flow chart illustrating the 
use of the simulation model for this scheme is shown in Fig. 12. 

With this technique, a forming tool radius is produced for each segment. A composite forming tool 
contour for the entire cross-sectional element is then determined by tangentially attaching the tool radii of 
adjoining segments In respective order, each having the length of the original segment. See Fig. 13. This 

45 creates the tool contour 220 needed to produce the cross-sectional element shape. To determine the 
surface of the three-dimensional forming tool surface, each forming tool curve is relocated in the plane for 
example, planes 205, 207, 209, from which the original cross-sectional element came. Once alt of the 
forming tool curves are in place, they are smoothly connected by fairing in a finished surface 221 (Rg. 14) 
between them with the aid of an electronic three-dimensional graphic system such as that known by the 

50 trademaric "CATIA". This finished surface 221 represents the definition of the forming tool contour. 

A second scheme exists for determining forming tool contours utilizing the age forming simulation 
model. Similar to the first scheme, again viewing Fig. 10B, the desired metal member 203 is analyzed with 
conveniently spaced apart planes 205, 207, 209 perpendicular to the tool reference plane 201 (Fig. 10A). 
Once again these planes are oriented across the contour and create cross sectional elements. Unlike the 

55 first scheme, however, the entire cross sectional element 211 is considered at once by the computer 
simulation. The element is first analyzed by considering its thickness profile. This can be done by 
considering the cross sectional element as a series of segments as in Fig. 11B. From this, a flat panel 
thickness profile 223 can be determined as seen in Rg. 15. This allows the panel to be described by a 
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coordinate system where the y-axis (ordinate) represents thickness (hi.h2,h3,h4, respectively) and the x-axis 
(abscissa) represents distance across the panel (h. b. b, U, respectively). This information is provided to 
the simulation model. Additionally, the shape of the cross sectional element 21 1 is described in terms of an 
x-y coordinate system. Also, the initial forming tool contour is approximated as having the same basic 
5 shape as the cross sectional element but containing smaller radii of curvature. The forming tool contour Is 
also described in terms of an x-y coordinate system. Both it and the cross sectional element coordinates 
are supplied to the simulation model. The measurement of the x-y coordinates of each of these shapes Is 
facilitated by the use of the electronic graphics package noted above. Alternatively, the shape of the cross 
sectional element and the shape of the forming tool contour could be described as a mathematical 
10 expression and be supplied to the computer simulation in that form. 

With the thickness profile and the shape description available, the computer simulation model is ready 
to operate. The model creates a simulated metal member of uniform width containing the thickness profile 
of the cross sectional element. The model simulates this created member being age formed against the 
approximated forming tool shape. The simulation concludes by describing the final shape of the created 
15 member after it experiences spring back. The final shape of the created member Is compared to the shape 
of the cross sectional element. If the initial approximation of the forming tool shape is correct, then no 
difference will be found between the shape of the simulated age formed member 225. represented by dash 
lines in Fig. 16, and the shape of the cross sectional element 211, represented by solid lines in Fig. 16. 
A comparison between these two shapes can be made by taking the vertical difference between them 
20 at various x-coordinate locations (Ai.A2,A3). also as shown in Rg. 16. Each difference can be compared to 
a predetermined tolerance value to establish whether or not the differences are acceptable. If the 
differences are not acceptable, then the simulated forming tool shape is adjusted and the age forming 
simulation is performed again using the same created member as before. The differences indicate the 
amount of underform or overform that occurs in the simulated member as compared to the cross sectional 
25 element. The simulated forming tool shape is adjusted by adding or subtracting the differences measured 
between the simulated age formed member and the cross section element. For instance, if the comparison 
indicates the simulated age formed member is underformed at a particular x-location. then the simulated 
forming tool is redefined as having more deflection at that particular location, the amount being equal to the 
difference between the simulated member and the actual cross sectional element. 
30 This iterative process continues until the differences measured are within the acceptable tolerance level. 
The simulated forming tool contour used to produce the acceptable created member shape is outputted 
from the simulation program in either an x-y coordinate format or as a mathematical expression. This 
represents the tool contour needed to produce the cross sectional element. 

With this information for each cross sectional element, the three-dimensional tooling is created in the 
35 same manner as described with respect to the first scheme. 

A third scheme for determining forming tool contour is very similar to the first scheme in the way it 
analyzes the complex shaped member to be formed. 

Conveniently spaced planes are utilized and cross sectional elements are created. Segments of each 
element are determined from sections having substantially similar radii and substantially similar thicknesses, 
40 each with a corresponding length. The great difference between the two schemes is in the manner in which 
the age forming simulation is utilized. In the first scheme, the simulation model is used to determine unique 
forming tool radii for each segment. In this scheme, the simulation model is used to prepare a general 
relationship between forming tool radii and formed members. This method is attractive in that simulations 
can be run for a given material prior to any knowledge of the shape of the complex member requiring form. 
45 The relationship is developed by simulating the forming of a range of thicknesses against a variety of 
forming tool radii where Rt is tool radius. Rp is part radius, and h is part thickness. This relationship can be 
displayed graphically as shown in Fig. 17, With this relationship, each segment can be analyzed according 
to its thickness and final formed radius. The resulting tool radii can be used to create a composite tool 
curve 227 for the entire cross sectional element similar to the first scheme shown in Fig. 13. The forming 
50 tool curves can be located in their respective planes, then faired together Into a surface to define the 
complete forming tool surface. 

Claims 

55 1- A method of forming a desired metal member having a surface contour of complex shape from an 
unformed member comprising the steps of: 

a) overforming the unformed member in a forming tool having a contour of smaller curvature than 
the contour of the desired member; 
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b) constraining the unformed member in the overformed condition; 

c) applying a thermal aging cycle to the member; 

d) cooling the constrained member following the thermal aging cycle; 

e) releasing the constrained member from the condition imparted by step (b) and allowing It to 
spring back to a dimensionally stable condition which defines the desired member having a surface 
contour of complex shape; 

(f) simulating the application of steps (a) through (e) to a geometrical representation of a member 
having the materia) properties of the desired metal member; and 

(g) determining from step (f) the contour of smaller curvature in step (a) to produce the desired 
metal member following step (e). 

A method as set forth in Claim 1 wherein step (f) includes the step of: 

(h) predicting the behavior of the desired metal member during steps (a) through (e) in terms of 
physical phenomena exhibited by the member including elastic deformation, Inelastic deformation, 
creep, and stress relaxation. 

A method according to claim 2 wherein step (h) includes the steps of: 

0) raising the temperature of the metal member according to the thermal aging cycle; 

(k) applying substantially uniformly distributed pressure on the metal member to conform to the 

contour of the tool resulting in elastic and Inelastic deformation of the metal member; and 

(I) holding the metal member against the tool for a predetermined time according to the thermal 

aging cycle whereby stress relaxation in the metal member occurs at a substantially constant strain; 

and 

(m) removing the uniformly distributed pressure from the metal member thereby allowing the 
remaining stresses in the member to decrease substantially to zero. 

A method according to claim 3 including the step of: 

(n) predicting the radius of curvature of the metal member resulting from both metallurgical and 
mechanical stress relaxation under deformation to achieve in the metal member the dimensionally 
stable condition which defines the desired metal member having a surface contour of complex 
shape. 

A mettiod according to any of claims 2 to 4 wherein in step (h) the physical phenomena exhibited by 
tiie desired metal member are defined by the following relationships: 

a + (E/T)(sinh(02/(E2f)))"E€ 

f = (a/r)((sinh(a2/(E2f)))". (sinh(/32f))n) 

f(0) = fo 

wherein: 

o = stress 
€ = strain 

f = internal variable 
o. iS. fo, n = material constants 
E = Young's modulus of elasticity 
r = relaxation time. 

A method according to any of claims 2 to 4 wherein in step (h) the mechanical and geometrical 
properties of the desired metal member are combined and defined by the following relationships: 
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6i = - ((EA)(sinh(o?/(^fi)))" + 3lV(2hS) 



-{24E/(th^) (slnh(a^/OE? mfyOyiiy/h)); l « 1,2 N 



ti = (a/T)((slnh(<^/(^fO))" - (sinh(P? ©)"); i = IX...^ 



< = (3?p/(2Eh^ - {24/fth^) I (slnhCo^/tg^OM^W/h 




wherein: 



a = Stress 

f = internal variable 

E = Young's modulus of elasticity 

T = relaxation tinne 

p = pressure 

h = thickness 

I = length 

y = y axis vertically upward (the undeformed neutral axis corresponds to y = 0) 
X = curvature 

/9, a. n = material constants. 

A method according to any of claims 1 to 6 wherein in step (f) the mechanical and geometrical 
properties of the desired metal member are combined and defined by the following relationships: 



a = Stress 

f = internal variable 

E = Young's modulus of elasticity 

T = relaxation time 

p = pressure 

h = thickness 

I = length 

y - y axis vertically upward (the undeformed neutral axis corresponds to y = 0) 
/8, a, n = material constants. 

A method according to any preceding claim including the steps of: 

(t) predicting a relationship between the radius of curvature of the contour of the desired metal 
member and the radius of curvature of the forming tool; 

(u) knowing the radius of curvature of the desired metal member, determining from the relationship 
of step (t) the radius of curvature of the forming tool. 

A method according to claim 1 including the steps of: 

(h) providing a geometrical representation of the desired complex shaped metal member; 




f| = (aAt)((sinh(a?/(^fO))" - (^(f fOf)l i » 1 A^-N 




wherein: 
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(i) passing a plurality of spaced apart planes through the geometrical representation of the desired 

member at spaced apart locations to thereby form a plurality of cross sectional elements: 

0) dividing each of the cross sectional elements into a plurality of segments, each segment having a 

length and a substantially uniform thickness and a substantially uniform radius of curvature; 

(k) determining from the simulation of step (f) a radius of curvature of the fomiing tool for each 

radius of curvature sought for each segment; and 

(I) from the tool radii calculated in step (k), developing tool curves for each cross sectioned element 
for each of the planes of step (i); 

(m) locating each tool cun/e determined in step (I) in each respective plane of step (i) in which its 
associated cross sectional element is located; and 

(n) joining the adjacent tool curves for all of the planes to thereby develop a surface contour for the 
tool. 

A method according to claim 1 including steps of: 

(h) providing a three-dimensional computer generated graphical representation of the desired 
complex shaped metal member; 

(i) passing a plurality of mathematically defined spaced apart planes through the graphical repre- 
sentation of the desired member at spaced apart locations to thereby form a plurality of mathemat* 
Ically defined cross sectional elements; 

(j) dividing each of the mathematically defined cross sectional elements into a plurality of mathemat- 
ically defined segments, each having a substantially uniform thickness and a substantially uniform 
radius of curvature; 

(k) determining from the simulation of step (f) a radius of curvature of the forming tool for each 
radius of curvature sought for each mathematically defined segment; and 

(I) from ttie tool radii calculated in step (k), developing tool curves for each of the mathematically 
defined planes of step (i) and thereby developing a mathematically defined surface contour for the 
tool; 

(m) locating each tool curve developed in step (I) in each respective plane of step (i) in which its 
associated mathematically defined cross sectional element is located; and 

(n) joining the adjacent tool curves for all of the planes to thereby develop a mathematically defined 
surface contour for the tool. 
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FIG. 4 
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FIG. 7 
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FIG. 8 
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FIG. 12 
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CONFORM TO FORMING 
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